SKYRME-RPA DESCRIPTION OF DIPOLE GIANT RESONANCE IN HEAVY 

AND SUPERHEAVY NUCLEI 
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The E1(T=1) isovector dipole giant resonance (GDR) in heavy and super-heavy deformed nuclei is 
analyzed over a sample of 18 rare-earth nuclei, 4 actinides and three chains of super-heavy elements 
(Z=102, 114 and 120). Basis of the description is self-consistent separable RPA (SRPA) using the 
Skyrme force SLy6. The self-consistent model well reproduces the experimental data (energies and 
widths) in the rare-earth and actinide region. The trend of the resonance peak energies follows the 
estimates from collective models, showing a bias to the volume mode for the rare-earths isotopes 
and a mix of volume and surface modes for actinides and super-heavy elements. The widths of the 
GDR are mainly determined by the Landau fragmentation which in turn is found to be strongly 
influenced by deformation. A deformation splitting of the GDR can contribute about one third 
to the width and about 1 MeV further broadening can be associated to mechanism beyond the 
mean-field description (escape, coupling with complex configurations). 



I. INTRODUCTION 

The isovector giant dipole resonance (GDR) is a most 
prominent and much studied excitation mode of nuclei, 
see e.g. [HQ. Nonetheless, it remains a subject of actual 
interest as there are many aspects which deserve more de- 
tailed investigations as, e.g., photo-excitation cross sec- 
tions in exotic nuclei which play a role in astrophysical 
reaction chains [3( or isotopic trends of the GDR includ- 
ing the regimes of deformed nuclei. The present paper 
aims at a theoretical survey of the GDR in isotopic chains 
of heavy and super-heavy nuclei. 

The high importance of the GDR has triggered since 
long many theoretical surveys analyzing the intriguing 
aspects of nuclear collective motion, starting from a 
purely collective description [1, [j| and slowly establish- 
ing a link to microscopic models in the framework of the 
Random-Phase Approximation (RPA) 0, 0]. The the- 
oretical description has much developed over the years. 
The majority of RPA investigations in the past employed 
shell model potentials plus an effective residual interac- 
tion (Migdal theory) @, ©, [rJ. In the meantime, self- 
consistent nuclear models have been steadily improving 
towards a reliable description of nuclear structure and ex- 
citations, for reviews see, e.g., [Ill [H, [lj| . These models 
belong to the class of density functional methods which 
aim at a universal energy functional for a given system 
from which all static and dynamics equations are derived 
in a strictly variational frame [lil . Il5j . Such density func- 
tional models, being rather universal by construction, are 
promising for investigation in exotic areas, e.g. for r- 
process, drip-line and super-heavy nuclei. The studies in 
this paper are based on the Skyrme functional which has 



been introduced in [l|| [l~7l a nd extended to a dynamical 
description shortly after 

im nana. 

The performance of 
self-consistent RPA calculations using the Skyrme func- 
tional had been tested systematically in [U [12] and it 
was found that one can have a reliable description when 
taking care to chose an appropriate parametrization. 

Systematic scans through the isotopic landscape and 
the study of exotic nuclei run over many deformed nu- 
clei. Fully fled ged RPA calculations for deformed nu- 
clei are feasible [23| , but extremely time consuming, not 
suited for systematic investigations. Super-heavy nuclei 
are especially demanding due to a coexistence of two ob- 
stacles, large size and deformation. Accurate but less 
demanding RPA techniques are needed. To that end, 
the separable RPA (SRPA) based on the Skyrme func- 
tional was recently developed [24], [2f| . The self-consistent 
factorization of the residual interaction in SRPA reduces 
the computational expense dramatically and so gives way 
to systematic explorations of nuclear giant resonances in 
both spherical and deformed nuclei [24|-[28|]. In this pa- 
per we concentrate on the isovector (T=l) electric GDR. 

As was shown in our previous studies [13|"[11|, SRPA 
provides an accurate description of the GDR in spheri- 
cal and deformed nuclei. We obtained good agreement 
with experiment for 154 Sm, 238 U and Nd isotopes with 
A=142, 144, 146, 148, 150. Eight different Skyrme 
forces were checked in these investigations. In this pa- 
per, we aim at a large systematics over the isotopic 
landscape and decide for one parameterization, namely 
SLy6 [29| which was found to provide a satisfying de- 
scription of the GDR for spherical and deformed nu- 
clei. In a first step, SRPA results for the GDR will 
be compared with all available experiment data in rare- 
earth and actinide regions. In particular, we consider 
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splitting) and their trends with system size will be ana- 
lyzed. 

In a second step, we will investigate the GDR in super- 
heavy nuclei for the three isotopic chains: Nobelium with 
Z=102 (A=242, 248, 254, 262, 270), Z=114 (A=264, 274, 
284, 294, 304) and Z=120 (A=280, 288, 294, 304, 312). 
As discussed below, this set covers most of the important 
mass regions and so is sufficiently representative. The 
main features of the GDR are analyzed and compared 
with those in rare-earth and actinide region. 

The paper is organized as follows. The calculation 
scheme, methods of analysis and choice of the isotopes 
are sketched in Sec. [TTJ In Sec. IHII the results of the cal- 
culations for the GDR in rare-earth, actinide and super- 
heavy nuclei are discussed. A summary is done in Sec. 

El ' 



II. FRAMEWORK 

The calculations are performed in the framework of 
SRPA [H| with the Skyrme force SLy6 [H which pro- 
vides a satis fying de scription of the GDR for heavy nu- 
clei [111, [5^, I27I l28|. The contribution to the residual 
interaction from the time-odd current density, Coulomb 
and pairing (at the BCS level) are taken into account 
[25l | . The calculations employ a cylindrical coordinate- 
space grid with the mesh size 0.7 fm. The calculation 
box has 24-35 mesh points depending on the nuclear size 
and deformation. For generators of the separable inter- 
action, we use four input operators, rYlfj, and j((iir)Yi^, 
i = 1,2, 3, with qi chosen following the prescription [24| . 
These operators are shown to give accurate results in 
spherical nuclei (24j. Besides, the operator r 3 Y 3fi is added 
to take into account the multipole mixing of excitations 
with the same projection p and space parity it. 

For all nuclei, the equilibrium quadrupole deformations 
are found by minimization of the total energy. The de- 
formations are characterized by the charge quadrupole 
moments Q2 and related dimensionless parameters 02 as 



dfp p (r)r 2 Y 20 , 
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(1) 



where p p (r) is the proton density in the ground state, 
< r 2 > p — J dfp p (r)r 2 / Z is the r.m.s. proton radius and 
Z is the number of protons. 

The calculations use a large basis space of single- 
particle states to expand the SRPA operators with two- 
quasiparticle states extending up to ~ 65 MeV, see Table 
I. The energy-weighted sum rule 
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is then exhausted by 92-95%. This sum rule includes 
the isovector effective mass m\ as the velocity-dependent 
terms are involved to the calculations, see discussion in 
[H. In SLy6 we have actually mf/m=0.80. 



TABLE I: Characteristics of the configuration space used in 
the calculations: minimal E m in and maximal E mam single- 
particle energies, number of the single-particle levels K for 
protons and neutrons, number of two-quasiparticle dipole 
states N2 qp (for branches p — and 1 altogether). See text 
for more details. 
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Since SRPA includes the dipole time-odd momentum- 
like operators [25| . the center-of-mass mode should, in 
principle, be placed correctly at zero energy. However, 
the finite computational box breaks translational invari- 
ance such that the spurious isoscalar El strength be- 
comes concentrated at the region of 2-3 MeV. That is 
safely below the GDR which we are studying here. 

The energy- weighted isovector dipole strength function 
is computed as 
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smoothed by the Lorentz function 
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with the averaging parameter A=2 MeV in most of the 
calculations. That averaging width was found to be opti- 
mal for the comparison with experiment and simulation 
of broadening effects beyond SRPA (escape widths, cou- 
pling with complex configurations). Further, 'J'o is the 
ground state, v runs over the RPA spectrum with eigen- 
frequency lo v and eigenstate >. 

To estimate the resonance energy centroid E, width T, 
and deformation splitting AE, the following prescriptions 
are applied. To determine E, the energy interval around 
the resonance with the strength larger than 10% of the 
maximal value is used and the centroid of the strength 
inside this interval is determined. The same method is 
implemented to find centroids Eq and E\ of p = and 
1 branches separately. Then the deformation splitting 
AE = Ei — Eq is obtained. The width T is determined 
at a half-maximum of the resonance. 

The experimental data for the GDR 0, HH, H [H, 
l35l . [H, [33, [H, Hi| include photoabsorption 

a = (7,n) + (7,p) + (7,np) + (7,2n) + (7,d) + ... + (7,/) , 

neutron yield 

a = (7, n) + (7, np) + 2( 7 , 2n) + 3( 7 , 3n) + ... + (7, /) , 
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FIG. 1: Isovector dipole strength in 156 - 160 Gd and 166 > 168 Er. 
The calculated strength (solid curve) is com par ed with the ex- 
perimental data for total photoabsorption [3ll . [32] (triangles) 
and neutron product [33j (rhombus). The branches of the res- 
onance with jj, = (left small) and n = 1 (right large) are rep- 
resented by dotted curves. The deformations are /?2=0.347, 
0.359 and 0.348, 0.353, respectively. 



and neutron product 

cr = (7, n) + (7, np) + (7, 2n) + (7, 3n) + ... + (7, /) . 

The photoabsorption data are preferable as involving all 
the main decay channels and so most corresponding to 
the strength function If the photoabsorption data 
are absent, the neutron yield and neutron product can 
be also used for a rough comparison since these data in- 
clude most of the main channels. However, one should 
take into account that both neutron yield and product 
omit (7, p) channel and so can underestimate the strength 
and change the resonance gross-structure. Besides, 
the neutron yield amplifies the neutron contributions 
(7, 2n), (7, 3n), ... and hence the right wing of the reso- 
nance. In what follows, we use experimental data for pho- 
toabsorption in 156 Gd [3l|, 160 Gd 168 Er El, 174 Yb 
El, 178 ' 180 Hf [H, 182 ' 1M ' 186 W 1 232 Th 1 
neutron yield in i70,i72,i76 Y b [3J 
neutron products (7, n) + (7, np) + 
17 fHf H and (7, n) + (7, np) + (7, 2n) + (7, /) in 
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In rare-earth and actinide regions, we consider all nu- 
clei for which reasonable GDR experimental data exist 
(for exception of Nd and Sm isotopes already explored in 
our previous papers [25j, |26j, |27|, |28| ) . In super- heavy nu- 
clei, we look at three isotopic chains: Nobelium 2=102 
(A=242, 248, 254, 262, 270), 2=114 (A=264, 274, 284, 
294, 304) and 2=120 (A=280, 288, 294, 304, 312). As 
can be seen from (46|, these chains cover most interest- 
ing mass and deformation regions. Indeed, they involve 
the onset (2=102), the center (2=114) and the upper 
end (2=120) of the super-heavy region. Every chain ex- 
tends through the whole neutron interval at a given 2. 



FIG. 2: Same as in Fig. 1 for 17 o,i72,i74,i76 Y b. The exper- 
imental data are for neutron yield [34| (closed circles) and 
total photoabsorption [3l[ (triangles). The deformation pa- 
rameters are /?2=0.350, 0.347, 0.340, 0.327, respectively. 



Different deformation regions are covered. For 2=102 
we deal with well deformed nuclei and small variation of 
the quadrupole deformation. The chains 2=114 and 120 
show strong variations of the deformation with a strong 
decrease with increasing TV, i.e. when moving towards 
the magic neutron number A^=184 [40] . The proton num- 
bers 2=114 and 120 are tentatively magic [H S3 such 
that neutron shell structure acquires a decisive weight for 
sphericity or deformation. 



III. RESULTS AND DISCUSSION 
A. Rare-earth and actinide nuclei 

Results of SRPA calculations for rare-earth and ac- 
tinide nuclei are presented in Figs. [UH Note that for 
reasons of comparison the calculated strength function 
is rescaled so as to correspond roughly to the maximal 
magnitude of the experimental cross section. Moreover, 
because of the insufficient accuracy of the model and ex- 
perimental resolution (see discussion in the previous sec- 
tion) we skip here the analysis of the fine structure which 
manifests itself mainly at the top of the resonance. The 
main attention will be paid the resonance energy cen- 
troids and widths. 

Figures [1] and [2] show an excellent agreement with ex- 
periment for Gd, Er and Yb isotopes. The agreement is 
less perfect for Hf, W and Os shown in Figures [3] and 
21 the calculated strength exhibits a slight (~ 0.5 MeV) 
down-shift in Hf and W and a larger (~ 1 MeV) up-shift 
in Os. It is known that Os isotopes are soft to oblate 
quadrupole deformation, which is confirmed by our cal- 
culations. However, as we checked, SRPA calculations on 
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FIG. 3: Same as in Fig. 1 for and 182 > 184 . 186 W. 

The experimental data are for total photoabsorption [3l| ] (tri- 
angles) and neutron product [H| (closed circles). The de- 
formation parameters for "MTs.isogf and 182,184,186 w are 

#2=0.330, 0.296, 0.287 and 0.260, 0.252, 0.247, respectively. 



FIG. 5: Same as in Fig. 1 for 232 Th and 234 > 236 ' 238 IJ. The ex- 
perimental data are for total photoabsorption [37J (triangles) 
and neutron product [381.I39T] (rhombus). The deformation pa- 
rameters for 232 Th and 234 > 236 . 238 U are f3 2 =0.256 and 0.279, 
0.286, 0.287, respectively. 
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FIG. 4: Same as in Fig. 1 for i 8 6,188,190,192q s ^ The experi _ 
mental data are for neutron yield [36] (rhombus). The defor- 
mation parameters are /32=0.222, 0.217, 0.195, 0.172, respec- 
tively. 

top of the oblate isomer do not improve agreement with 
the experiment. The discrepancies for Os may be partly 
caused by using the neutron yield experimental data. As 
was discussed above, the neutron yield alone can am- 
plify the right GDR flank, thus resulting in some appar- 
ent up-shift as compared with the total photoabsorption 
cross section. Results for the actinides in Fig. [5] also 
look encouraging. Modest deviations in gross-structure 

of 234,236TJ 

can be explained by using the neutron product 
experimental data. 

It is worth noting that the comprehensive analysis of 



experimental data reveals noticeable (sometimes signifi- 
cant) deviations in GDR measurements of different ex- 
perimental groups (42| . Taking into account these uncer- 
tainties in the data, one may consider for all cases above 
the agreement with experimental data as very satisfying, 
at least for energy centroids and widths. 

The correlation of the quadrupole moments JTJ) with 
some resonance characteristics (width T, deformation 
splitting AE, energy centroid E) and trends of these 
characteristics with the mass number A are shown in Fig. 
|51 All nuclei in the sample have a significant quadrupole 
deformation. The calculated quadrupole moments are 
in excellent agreement with the experiment (43J. In the 
rare-earth nuclei the moments have a maximum in the 
middle of the region. Note that the dimensionless defor- 
mation parameters P2 as indicated in the captions of the 
previous figures are maximal at the onset of the region. 
The difference between Qi and maxima is related with 
the nuclear mass factor in . The direct contribution of 
the deformation splitting AE to the resonance width is 
maximal (~ 40%) in the first half of the region (A < 176) 
and then slowly decreases to 37 — 34% for Hf, 34 — 31% 
for W and 31 — 24% for Os. Furthermore, AE is a bit 
increased in actinides where it reaches 30 — 33%. This 
trend obviously correlates with /?2- See also the detailed 
discussion of the GDR width and structure in section 

Unci 

In Figure [5J:, the resonance energy is compared with 
empirical estimates based on Steinwedel- Jensen (SJ) [f| 

E SJ = 81A~ 1/3 MeV (5) 

and Bcrman-Fultz (BF) 0, Q 

E BF = (31.2A~ 1/3 + 20.6^ 1/6 )MeV (6) 
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FIG. 6: Calculated characteristics of rare-earth and actinide 
nuclei as function of mass number. Upper panel: Quadrupole 
moments Qi compared with experimental values (crosses) 
[43l ]. Middle panel: Widths V (upper set) and deformation 
splittings AE (lower set). Lower panel: Energy centroids as 
compared with the estimates $5$ (dotted curve) and ([6]) (solid 
curve) . 
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collective models. BF takes into account both volume 
and surface contributions and treats the dipole resonance 
as a combination of Steinwedel- Jensen Q and Goldhaber- 
Teller (GT) [i| scenarios. The calculated energies are 
closer to the SJ estimate ([5]) in the rare-earth region and 
to the BF estimate in actinides. As shown below, 
the BF estimate is also more appropriate in super-heavy 
nuclei. This agrees with a commonly accepted view that 
in heavy nuclei neither the density gradient (~ A -1 / 3 ) 
nor the nuclear surface impact (~ A^ 1 / 6 ) dominate the 
restoring force [2, 13, Hj| . 

B. Super- heavy nuclei 

The agreement of SRPA results with the experimental 
data in rare-earth and actinide nuclei encourages its fur- 
ther application to super-heavy nuclei. SRPA results for 
super-heavy nuclei are exhibited in Figs. ITlfTTJl 

Fig. [7] indicates that the GDR in this region is gen- 
erally similar to that in rare-earth and actinide nuclei. 
In particular, the resonance width correlates with the 
quadrupole parameter fa. Furthermore, the middle pan- 
els in Figs. [8]fT0l show that the direct contribution of 
the deformation splitting AE to the resonance width T 
does not exceed 40%, as in rare-earth and actinide nu- 
clei. Note that our quadrupole moments Q2 from the 



FIG. 7: Same as in Fig. 1 for isotopes of super-heavy 
nuclei No and Z=114 and 120. The dimensionless proton 
quadrupole deformation is indicated in the plots. 



self-consistent calculations agree nicely with the values 
obtained within the macroscopic-microscopic model [4(| , 
see the upper panels in Figs. [8lfTfJI The agreement persist 
even in mass regions with large variations of deformation. 

At the same time, the GDR in super-heavy nuclei 
shows some peculiarities. First, unlike the rare-earth nu- 
clei, its energy is much closer to the BF estimate ([6]) 
than to ([5]), which supports once again the treatment of 
the GDR in heavy nuclei as a mixture of SJ and GT 
modes. Maximal deviations from both estimates emerge 
at the ends of the isotopic chains, which is natural since 
these estimates do not parameterize any isospin depen- 
dence. Second, as is seen from Figs. [51ITU1 the decrease 
of the resonance energies with neutron number N levels 
off and is even reversed to some increase at the end of 
every isotope chain. This can be explained by increase 
of the symmetry energy E cx E sym = a sym {N — Z) 2 /A 2 
at these neutron-rich edges. Note that such a turnover is 
absent for lighter nuclei, e.g. in the chain of Nd isotopes 
(A=134-158) explored earlier in [28|. Probably this is 
because Nd isotopes do not reach so large neutron excess 
as the super-heavy elements. 
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FIG. 8: Calculated characteristics of No isotopes as a function 
of their neutron number. Upper panel: Quadrupole moments 
Qi (black squares) as compared with the values [46l ] (open 
circles). Middle panel: Widths Y (black squares) and de- 
formation splittings AE (black circles). Lower panel: Energy 
centroids (black squares) as compared with the estimates Esj 
§5§ (dotted curve) and Ebf © (dotted curve with crosses). 
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FIG. 9: Same as in Fig. 8 for Z=11A isotopes. 
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FIG. 11: The isovector dipole strength calculated with dif- 
ferent averaging in the strength function <(3j : A =0.5 MeV 
(dash curve), 1 MeV (thin solid curve), and 2 MeV (bold 
solid curve) . The parameters of the proton quadrupole defor- 
mation are given for every nucleus. 



C. GDR width and structure: general discussion 

In Figure fTTl the isovector dipole strengths calculated 
with different Lorentz averaging parameter A are com- 
pared for a representative set of nuclei. It is seen that 
smaller averaging, A =0.5 and 1 MeV, yields more fine 
structures, mainly at the resonance peak, and leaves the 
total width almost unchanged. As a rule, rJT = T(A = 
2 MeV)-r(A = 0.5 MeV) < 1 MeV. Since T(A =2 MeV) 
well reproduces the experimental widths, one may asso- 
ciate the difference ST to the smoothing effects omitted 
in the present RPA calculations (coupling with complex 
configurations, escape widths). In fact, the averaging 



with A —2 MeV effectively mimics these effects. Since 
ST <c T(A = 2AleV), the deformation splitting and Lan- 
dau fragmentation (distribution of the collective strength 
between nearby two-quasiparticle states) obviously dom- 
inate the total width. We estimate their contribution by 
70-90%, depending on the nucleus and its shape. 

Figure [TT] also shows that the averaging A —2 MeV 
chosen in our calculations is indeed most suitable for the 
comparison with the GDR experimental data (at least 
for the heavy nuclei considered here). This averaging 
does not cause significant artificial increase of the reso- 
nance width and, at the same time, allows to suppress 
the structure details which, in any case, are not resolved 
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TABLE II: Calculated RPA widths F and Fi of the resonance 
branches /i = and 1, the sum r + and the total width F. 
For every nucleus the deformation parameters 02 are given. 
The averaging is A = 2MeV . For more details see the text. 



Nucleus 


P2 


Widths [MeV 








To 


Ti 


To +Ti 
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156 Gd 


0.347 


3,05 


4,74 


7,79 


7,69 


i72 Yb 


0.347 


2,54 


5,08 


7,62 


7,28 


186 Os 


0.222 


2,65 


5,09 


7,74 


6,39 


238 tj 


0.287 


2,62 


5,11 


7,73 


6,86 


254 No 


0.308 


2,44 
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FIG. 12: The isovector dipole strength calculated with 
smoothing A =2 MeV in deformed 156 Gd and spherical 
304 120. The upper panels show full SRPA results (with resid- 
ual interaction) and the lower panels pure two-quasiparticle 
(2qp) strengths (without residual interaction). The separate 
branches /x=0 (small) and l (large) are plotted by the dot- 
ted line and their sum by the full line. 



in the experimental strength distributions. 

Besides the splitting of the GDR into two branches, 
the deformation also results in a considerable redistribu- 
tion of the strength within every branch. In other words, 
the deformation severely influences the Landau fragmen- 
tation itself. This effect is illustrated in Table II and 
Fig. [T2j Table II provides the widths of the resonance 
branches /i = and /i = 1 for a selection of nuclei. There 
is a large difference between spherical and deformed nu- 
clei. In spherical nuclei ( 304 114, 304 120, 312 120), we have 
To ~ Ti « r, while in deformed nuclei, Ti > Tq and their 
sum To+ri roughly covers the total width T. So, we see a 
strong deformation effect within the branches themselves. 
In order to analyze it, we plot in the lower panels of Fig. 
[I2]the unperturbed two-quasiparticle (2qp) strengths. In 
deformed nuclei ( 156 Gd) the inequality Ti > Tq appears 
already in 2qp strengths, which shows that this is simply 
an effect of the density of states. The residual interac- 



tion (SRPA case) does not cause essential changes in the 
relation between /i = and \l = 1 branches and preserves 
Ti > To in deformed nuclei. Besides showing the influ- 
ence of deformation, Fig. [T2]also illustrates the collective 
shifts from the unperturbed strengths to the final ones. 
The size of the shifts is, of course, related to the actual 
force SLy6. The shifts are 4.0 and 3.4 MeV in 156 Gd and 
304 120, respectively. 

In most of the nuclei the calculations indicate a small 
shoulder at the right flank of the resonance. The heav- 
ier the nucleus, the stronger the shoulder. At a small 
averaging width, the shoulder becomes more pronounced 
and may even show up as a separate peak, see e.g. results 
for A=0.5 MeV in Fig. [TTJ This effect is often absent, 
at least much less pronounced, in the experimental data. 
The shoulder takes place in both deformed and spherical 
nuclei and so is independent of deformation. It persists 
not only for SLy6 but also for most of other Skyrme forces 
[25l [26j . As was shown in [26| , the effect is caused by spe- 
cific 2qp structures with high angular momentum (thus 
large statistical weight) lying at the right GDR flank. It 
is to be noted that RPA neglects some broadening mech- 
anisms (escape widths, coupling with complex, 2p-2h, 
etc, configurations) which could soften these structures. 
The shoulder can be further enhanced if the Skyrme force 
overestimates the dipole collective shift [25, 26] . The case 
calls for further detailed exploration. 



IV. CONCLUSIONS 

The isovector giant dipole resonance (GDR) is system- 
atically investigated in rare-earth, actinide and super- 
heavy regions. The study covers 37 nuclei altogether. 
Mainly axially deformed nuclei are considered. In all 
the nuclei, the calculated quadrupole moments correctly 
reproduce the experimental data (rare-earth and ac- 
tinide regions) [43| or macroscopic-microscopic estimates 
(super- heavy region) [4^]. The calculations are per- 
formed in the framework of the self-consistent separable 
RPA approach (SRPA) based on the Skyrme functional. 
The force SLy6 is used. 

A satisfying agreement of the SRPA results with the 
available GDR experimental data is found for 22 rare- 
earth and actinide nuclei. Resonance energies and widths 
are well described. The trends of the peak energies are 
compared with simple estimates from collective models. 
It is confirmed that the Steinwedel-Jensen (SJ) model 
performs well for medium heavy nuclei while a mix of SJ 
and the Goldhaber- Teller scenarios is more appropriate 
for heavy nuclei. Encouraged by these results, the survey 
is extended to super-heavy nuclei where GDR in isotope 
chains with Z=102, 114 and 120 are explored. The GDR 
in super-heavy nuclei is found to behave similar to that 
in rare-earth and actinide nuclei. The peak energies are, 
again, better described by the mixed collective model, 
continuing the trend from heavy nuclei. A new feature 
in the super-heavy region is that the peak energies turn 
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from decrease to increase towards the neutron rich ends 
of the isotopic chains (close to the drip lines). 

We also analyze the GDR widths T. They are shown 
to be strictly dominated (at least 70-90%) by Landau 
fragmentation and deformation contributions. The di- 
rect deformation contribution through the splitting of 
the GDR into /i = and 1 branches achieves 40%. The 
Landau fragmentation is severely affected by the defor- 
mation as well, which modifies the branch widths and 
lead to Ti > L . The final step to agreement with ex- 
perimental pattern is achieved by Lorentz averaging the 
SRPA results, thus simulating missing broadening mech- 
anisms (e.g., escape widths and coupling with complex 
configurations). A modest additional broadening of ~ 1 
MeV suffices to reach realistic pattern. 
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